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THE AXIS RATIO DISTRIBUTION OF X-RAY CLUSTERS OBSERVED BY XMM-NEWTON 

Hajime Kawahara 1,2 

ABSTRACT 

We derive the ax is ratio distribution of X-ray clusters using the XMM-Newton catalogue 
(|Snowden et al.l l2008f ) . By fitting the contour lines of the X-ray image by ellipses, we confirm the 
X-ray distribution is well approximated by the elliptic distribution with a constant axis ratio and 
direction. Wc construct a simple model describing the axis ratio of the X-ray gas assuming the hydro - 
static equilibrium embedded in the triaxial dark matter halo model proposed bv lJing fe Sutol (|2002f ) 
and the hydrostatic equilibrium. Wc find that the observed probability density function of the axis 
ratio is consistent with this model prediction. 

Subject headings: X-rays: galaxies: clusters, cosmology: observations, galaxies: clusters: general 



1. INTRODUCTION 

Galaxy clusters have been playing an important role 
in the determination of cosmological parameters such as 
the matter density, the amplitude of the initial fluctua- 
tion, the Hubble constant. Upcoming surveys of galaxy 
clusters via optical, X-ray, lensing and the Sunyaev- 
Zel'dovich effect will enable one to probe new aspects of 
cosmology such as t he nature of the dark energy (e.g. 
lAlbrecht etafl 120061 ) and the initial non-Gaussianity. 
Therefore, the physical modeling of galaxy clusters is im- 
portant to properly interpret the cluster data. In particu- 
lar, non-sphericity adds significant uncertain ty to various 
cosmological applicatio ns of the clusters (e.g. lOguri et al.l 
12001 lDsJld~etgj1l2T)0l . 

The statistical nature of the non-sphericity has been 
inve stigated using N-body simulations by many authors 
(e.g.lJing fe Sutoll2002l:lKasun fe EvraxcKOTjl iPaz et al.l 



120061: lAllgood et al l 120061: iBett et al.l I2007D . Among 
them. iJing fe Sutol ( 20021 hereafter JS02) proposed a tri- 
axial model of the dark matter halo based on the detailed 
analysis of N-body simulations. They also provided the 
fitting formula of the probability distribution of the axis 
ratio. 

Their phenomenological model has been applied to 
many studies that investigate the systematic errors of 
the cosmological applications with clusters. The obser- 
vational test for the axis rati o distr ibution, however, is 
quite limited. lEvans fe Bridle! ()2009| ) analyzed 4281 clus- 
ters in the Sloan Digital Sky Survey using the gravita- 
tional lensing. By stacking all the clusters, they were 
merely able to conclude that the axis ratio of the pro- 
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Due to large 



uncertainty, however, the measurement of the distribu- 
tion of the axis ratio was difficult. Future surveys will 
make a significant improvement of the direct measure- 
ment of dark matter h aloes with new techniques (e.g. 
lHawken k Bridle! 120091) . 

On the other hand, X-ray clusters also represent 
the non-sphericity, which makes the large statistical 
error to measurement of the Hubble constant (c.g 
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iKawahara et al.ll2008al ). The theoretical predictions of 
the probability density function (PDF) of the axis ratio 
using X-ra y clusters have been considered by sev eral au- 
thors (e.g. IWang kFa^&M ILee fc Sutdl200lh . These 
are based on the assumption that the non-sphcricity of 
X-ray clusters originated from the non-spherical gravita- 
tional potential due to the dark matter halo. It is cru- 
cial to see if the non-sphericity of X-ray clusters can be 
explained by the non-sphericity of the underlying dark 
matter halo. 

There are a lot of papers which measured the 
axis ratio of the X-ray clu s ters observed by E i nstein 
(e.g. iMcMillan et al.l 119891: iBuote fe Canizaresl 119921: 
Mohr et all I1995D. ROSAT (e.g. iBuote k Canizared 
19961: IWang fc Ulmerlfl997l: iKolokotronis et al.ll2001ft an d 
Chandra (e g Hashimoto et al.l 120071: iFang et al.ll2009D . 
iFlores et al.l (|2007h have constructed the theoretical pre- 
diction of the gas distribution and justified their model 
by comparing with 8 simulated clusters. Applying this 
model to 46 simulated dark matter haloes, they predicted 
the average and variance of the axis ratio of X-ray clus- 
ters. Comparing with ROSAT samples they have found 
good agreement between the predicted and observed av- 
erage and variance. 

Because the XMM-Newton has the largest effective 
area among present X-ray satellites and its field of view 
is also large, one can investigate the details of the X-ray 
shape or isocontour lines of clusters including the radial 
dependence of the axis ratio and the direction of semi 
major axis for each cluster. While the axis ratio con- 
straints are not expected to be much better than have 
been obtained previously, it is useful to have an indepen- 
dent computation to serve as a check on systematic errors 
between different data sets and different computational 
methods. In this paper, we perform detailed compar- 
isons of the PDF derived by the X-ray analysis using high 
quality data set of th e XMM-Newton cluster catalogue 
(|Snowden et al.ll2008l ) with the simple theoretical predic- 
tion of the PDF based on the phenomenological model 
of shape of underlying dark matter haloes provided by 
JS02. 

The rest of the paper is organized as follows. In §2, we 
describe the cluster data we use and the data process- 
ing such as fitting clusters and estimating the statistical 
errors. In §3, in order to model the non-sphericity, we in- 
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vestigate the radial dependence of the axis ratio and the 
axis direction. After describing a simple model of the gas 
non-sphericity, we compare the non-spherical properties 
of X-ray clusters with the theoretical prediction. Finally, 
we summarize our results in §4. 

2. METHODS 

We anal yze 70 clusters in th e XMM-Newton cluster 
catalogue iSnowden et al.ll2008l ). The surface brightness 
map of each cluster is normalized by the exposure map 
and has the energy range of 0.5-10.0 keV band. Al- 
though XMM-Newton has three CCD imagers, 2 MOSs 
and pn, we simply combine the images of MOS1 and 
MOS2 and do not use the data of pn. Based on the 
light curve with the energy range of 10-15 keV, we ex- 
clude the period of " flares " due to the soft proton 
background from our analysis. The point sources are 
also excluded by the wavelet detector in the CIAO pack- 
age. The images are binned in a pixel size of 1 % 
of the virial radius r2oo- The virial radius and mass 
are estimated by the relations found in the simulation, 
r 200 = 1.25 x (1 + z)" 3 / 2 (r/10.0kcV) 1 / 2 fe- 1 Mpc and 
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1.5 x (f/10.0key) 3 / 2 / t- 1 10 15 M^ (lEvrard et al.l 
119961: 1 Arnaud fc Evrardl 1 19991 ) . We adopt the averaged 
cluster temperature in the 1' — 4' annulus given by 
ISnowden et all (|20M ). 

We perfo rm the elliptical fitting according to 

iJedrzeiewskil |l987). The elliptical fitting has four free 
parameters, the minor to major axis ratio, q, the direc- 
tion of the long axis 0, and the center r c . In §??, we 
compare the shape with the prediction based on the re- 
sults of JS02 and they measure their PDF around the 
radius of ~ 0.3r2oo- Therefore we fix the semi-major axis 
R t of the ellipse to 0.1, 0.2, 0.3 and 0.4r 200 for i = 1,2,3 
and 4. 

The background level Sxb 1S estimated at the outer 
part of the surface brightness profile, Sx around the cen- 
ter. We compute the average of signal to noise ratio on 
the ellipse as 
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where N B(Ei) is the average count of the background 
estimated by Sxb an d the exposure map at Ei. we adopt 
the number of sam pling points, iV sam p = 360. We use the 
ellipse that satisfy S/N > 1. 

In order to sec the radial dependence of shape we inde- 
pendently fit the center as free parameter for each con- 
tour line). Because there is no cluster with strong dou- 
ble peaks fortunately, except for RXJ 658-55 (the Bullet 
cluster), the best fit centers are close to the maximum 
surface brightness point in almost clusters. The statisti- 
cal errors of q and due to the Poisson noise are esti- 
mated by the Monte Carlo method. We create 100 mock 
clusters by redistributing the photon number according 
to Poisson distribution. The standard deviation a q {Ri) is 
derived by fitting the contour lines of each mock sample 
with the same fitting method. 

Figure [1] displays an example of the elliptical fitting. 
We exclude parts of contour lines in the masked region 
from the fitting. There are 61, 56, 39, and 13 clus- 
ters that satisfy equation (TT]) for R/r 2 oo = 0.1,0.2,0.3 



and 0.4. We exclude 7 clusters, Coma, Perseus, A1060, 
M87, A2256, A3888, and RXJ0658-55 from our analysis 
because we cannot fit even -R/V200 = 0.1 due to their 
rcdshifts or the virial radius, their position in the field of 
view or the bullet cluster. The best fit value of q and its 
error for each cluster are listed in Table [T] The errors of 
q, a q is typically a few percent of q. 
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Figure 1. An example of the elliptical fitting of the X-ray clus- 
ter (A2597). Black curves are the fi tting ellipse f or R/ r2oo = 
0.1,0.2,0.3 and 0.4 by the method of IJedrzeiewskil (19871 ) ■ Four 
white dashed and dotted lines indicates the direction of the semi 
major and semi minor axes, respectively. 



3. PROBABILITY DISTRIBUTION OF THE AXIS RATIO 

Figure [5] shows the radial dependence of the projected 
axis ratio and the axis direction. In the upper panel, 
we plot the average of q for clusters as a function of 
R/f200- Although there is a slight increment of the axis 
ratio with z, the difference between R/r2oo = 0.1 and 
0.3 does not significantly change (w 0.04). The lower 
panel shows the radial dependence of the difference of 
9(i?i) around the radial averages, A0(i?j) = 0(i?i) — 
0. The average axis direction is computed from_thc 
average of the unit vector with an angle <d(Ri), = 
arccos [u x /{u x + m 2 ) 1 / 2 ], where u x = ^\cos0(-R,) and 
u y = £\ sin0(i?i). For most clusters, the axis directions 
are aligned within 10 degree. The points with A0 > 
7r/8 have large errors as displayed by color. These large 
statistical errors of the axis direction mainly originate 
from the uncertainty due to their high axis ratio (near 
spherical) . Despite some variations of the axis ratio and 
the axis direction, we find that the surface brightness 
contours are well approximated by the ellipses with the 
same direction and axis ratio for different radii. 

JS02 reported the negative redshift dependence and 
the (slightly) positive mass dependence on the axis ra- 
tio. Figure [3] displays the redshift and mass dependence 
of the axis ratio. Since the present XMM-ncwton cata- 
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Figure 2. Radial dependences of the projected axis-ratio and the 
axis direction. The upper panel shows the average profile of the 
axis ratio q. We average over 61, 56, 39, and 13 clusters for 
R/f200 — 0.1,0.2,0.3 and 0.4, respectively. Black errors indicate 
the standard deviation for clusters. The Lower panel display the 
difference of thc_axis direction, AO(Rj) = 0(i?;) — 0. The average 
axis direction is derived by the average of the unit vector with 
0(i?i), = arccos [u x /(u x + u 2 ,) 1 / 2 ], where u x = J2i cos @(Ri) 
and u y = J^. sin 0(Rj). Each dashed line connects A0 for different 
radii of an identical cluster. Each color indicates the standard 
deviation of 0(i?i) for the j-th cluster ffe,j : re d, orange, green, 
cyan, and blue. 

logue is heterogeneous, Figure [3] docs not show any clear 
dependence of the axis ratio on z or M 2 oo- 
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Figure 3. Scatter plot of clusters with different axis ratio on a 
Z-M200 plane. Different types of symbols indicate the axis ratio of 
the cluster: 0.9 < q (circle), 0.8 < q < 0.9 (cross), 0.7 < q < 0.8 
(triangle), and 0.6 < q < 0.7 (square). 

Let us compare the observed axis ratio distribution of 
X-ray halo with a theoretical expectation. We first con- 
struct a simple model of the axis ratio distribution of the 
X-ray halo based on the triaxial model of the dark mat- 



ter halo. Our assumptions are as follows, 1) The dark 
matter halo is a homcoid ellipsoid, (has the constant axis 
ratio and perfect alignment of the axis direction) . 2) The 
gas distribution is perfectly governed by the dark matter 
potential through the hydrostatic equilibrium. In other 
words, we neglect the self gravity of the gas and regard 
a cluster as a completely relaxed system. 3) The X-ray 
surface brightness distribution depends on the gas den- 
sity only, that is, we assume that the cluster is isother- 
mal. Although the actual simulated halo shows a slight 
radial dependence of the axis ratio as reported by JS02 
(see Figure 3 of their paper), we assume the constant 
axis model. 

JS02 suggested the triaxial density model with three 
axes a\ > a 2 > 03 described as 

P(R) = .... a , P " 7Z -SZ5- , (2) 



(R/R s f' a ] 



where a = 1 for the cluster size halo (the triaxial NFW 
profile), a = 1.5 for the galaxy size halo and R is the 
major axis length: 
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where \i A = 03/01 and [Ib = a 2 /a\ denote the 
minor-to-major and medium-to-major axis ratio, re- 
spectively. JS02 also provided the fitting formula of 
the joint probability density function (PDF) of [La 
and [i B at ^4 (3) = 2500, Pdm(/"a, fJ-B] M 200 ; z) = 
Pa{^a; M 200 ; z)pb(mb|^a), where 
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where M* is the characteristic nonlinear mass at z so 
that the rms top-hat smoothed over density at the scale 
is 1.68, a s = 0.113, and r m - ln = fiA for \xa > 0.5 and 
r m in = 0.5 for ^a < 0.5. The is R with which the 
ellipsoid having overdensity 2500. Note that A(3) = 2500 
corresponds to the typical radius « 0.3r2oo- 

Under the hydrostatic equilibrium assumption, the gas 
isodensity surface is identical to the isopotcntial surface. 
The potenti al of a homeoid at t he po sition [x, y, z) is 
provided bv lBinnev fc Tremainel (|1987t) . 
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Assuming the density distribution in equation @, we 
obtain 
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We can rewritten the potential as, 



$(x',y',z') = -2iTG 
Z(x',y',z<) 



R 2 p 



p A p B ^(x',y',z') 



/-[i + K)- 1 !"- 8 dT , {7) 

o V(l + r)(^+r)(^+r) 

(8) 

where x',y', z' and ml denot e x, z and ro norm alized by 
R s , respectively. Although ILee fc Sutol (|2003[) provided 
a perturbative approach to solve $, we directly derive 
the axis ratio by numerical integration. The isopotential 
surface defined by §(x, y, z) = const = <J>(:z.*, 0, 0) can be 
approximated by the triaxial ellipsoid with the minor to 
major and the medium to major axis ratio, 

(9) 
(10) 



fi>g,A(fi>A, Pb', = — r, and 



where y^ and z 1 ^ are obtained by solving following rela- 
tions, 

S«,0,0)=S(0,0,4) = const, (11) 
S(ai I 0,0) = S(0,i/.,0) = const. (12) 

Thus, we regard that the gas density distribution as an 
ellipsoid with constant axis ratios p g ,A and p g .s- 

One uncertainty of this model is the choice of . Be- 
cause the fitting formula given in equation is based 
on R ss 0.3 7-200, it is natural to adopt x'„ = 0.3 r2oa/R s = 
0.3 c, where c is the concentration parameter. For sim- 
plicity, we assume the typical value c = 3 found in simu- 
lations (JS02). Because most contour lines we investigate 
are below 0.3r2oo, we also consider x' t = 0.3,0.6 which 
approximately corresponds to the contours i = 1 , 2 in 
addition to 0.9 (i = 3). 

The axis ratio of the ellip se made by th e projection of 
the ellipsoid is obtained by (|Binnevlll985h 
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where 9 and <f> are the polar coordinates of the line of 
sight. The PDF of the projected axis ratio is given by 
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where 5d{x) is the delta function. As shown in Figure 
02 the cluster sample is not homogeneous for mass and 
redshift. We do not know the mass or redshift selection 
function of the catalogue. Therefore we use the ensemble 
averaged PDF, not using the fixed mass and redshift, 
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Figure 4. Comparison the theoretical prediction of the axis ratio 
distribution (lines) with observation (error bars). We use the ra- 
dial averaged axis ratio of each cluster. Error bars are Poissonian. 
Lines in the top and middle panels are theoretical curves with the 
density profile (Eq.[2]) with a = 1 and a = 1.5. Solid, dashed, 
and dot-dashed lines indicate x' t = 0.9, 0.6 and 0.3, respectively. 
For reference, we plot the axis ratio distribution of the underlying 
dark matter halo, which is expected from the lensing observation 
for instance. The bottom panel displays the PDF for different 
radii. Black, red, green, and blue lines indicate R = 0.1,0.2,0.3 
and 0.4r200, respectively. 

for comparison with observation. We adopt the virial 
mass and redshift of each cluster for M2oo,i and Zj. 

Figure |4] displays the observed axis ratio distribution 
and the theoretical prediction described in equation (|15|) . 
Because the radial dependence of the ellipse is relatively 
small, we adopt the radial average for each cluster as the 
observed axis ratio. We plot the theoretical prediction 
described in equation (|T5|) by performing Monte Carlo 
methods. The number of realization for each cluster is 
128 2 . The top and bottom panels show the results for 
a = 1.0 and 1.5, respectively. We find that both models 
significantly agree with observational result although the 
NFW triaxial model is slightly better than that of a — 
1.5. We also confirm that the model uncertainty due to 
the selection of x'^ is relatively small. The bottom panel 
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displays the PDFs for the different radii R = 0.1, 0.2, 0.3, 
and 0.4r2oo- There is no strong dependence of the PDF 
on the radius between 0.1 — 0.4r2oo 

Because our model is quite idealistic, many sources of 
the systematic errors are possible. For instance, the devi- 
ation from the isothermality affects the relation between 
the gravitational potential and X-ray surface brightness. 
Under the hydrostatic equilibrium assumption, the con- 
tour of the gravitational potential is identical to nT and 
n 2 A(T). Let us consider two points A and B on the isopo- 
tential surface. The relation of the local X-ray emissivity 
at A (e A ) and B (e B ) is e B = [T%A(T B )/(T%A(T A ))]e A « 
(Ta/Tb) 3 ' 2 £a, where T A and T B are temperature, n A 
and n B are gas density at A and B, and A(T) is the 
X-ray emissivity. If assuming the j3 model, the surface 
brightness profile follows oc r~ 3/3 at out of the core. Then 
we approximately estimate the change of gas emissiv- 
ity contour from as A - (T B /T A ) 3 / 4 . Therefore, 10 % 
global deviation of temperature from isopotencial surface 
corresponds to < 10 % change of the iso-emissivity sur- 
face or the axis ratio of X-ray haloes. We also note that 
our model use the shape of the potential at a given ra- 
dius and ignore the projection effect of the gas density. 
This assumption is also a possible source of systemat- 
ics. Moreover, the hydrostatic equilibrium assumption 
itself may not valid for merger clusters. A rich variety 
of structural comp lexity seen by X-ra y observations like 
shock fronts (e.g . . I Jones et al.l l200"2t ). col d fronts (e.g., 
Markevitch et al.l 12000^ 7 X- ray holes (e.g.. iFabian et al.1 
2002f ) . gas fluctuations (e.g. JKawahara et al.l2008bh . and 
AGN feedback are also possible sources of the systematic 
error of the theoretical prediction. Considering the sim- 
plicity of our current model, however, the agreement is 
remarkable. 

3.1. Comparison with Previous Studies 

In this section, we compare our results with other 
studie s by Einstein, ROSAT, and Chandra. iMohr et al.l 
(| 1995T) provided the distribution of the axis ratio for 
51 clusters observed by Einstein. The mean axis ra- 
tio of their data ((q) = 0.80) is similar to our results 

«<Z) = ESx(EfcT?,'(c«)/w)/JVd - 0.791 ± 0.004), 
where N c \ = 62 is the number of the clusters and i max 
is the number of contour lines for each cluster. While 
the dispersion of their data (s q = 0.12) is slightly larger 
than ours (s q = 0.095 ±0.007), it is likely that the spatial 
resolution and sensitivity between Einstein and XMM or 
the different fitting methods make a difference of the dis- 
p ersion. 

iFlores et al.l (I2007D discussed an average cllipticity (e = 
1 — q 2 ) and its dispersion by the flux- limited samples 
of ROSAT clusters. They provided the mean cllipticity 
(e) = 0.376 ±0.019 and its dispersion s e = 0.122 ±0.014. 
Although the mean cllipticity of our analysis ((e) = 1 — 
(<7 2 ) = 0.365 ±0.007) is slightly smaller than their result, 
the dispersion (s e = .144 ± 0.005) is good agreement. 
lHashimoto et al.l (|2007[ ) also derived the mean ellipticity 
e' = 0.22 with 101 clusters observed by Chandra (their 
ellipticity is defined by e' = 1 — q) . Converting to the 
mean axis ratio, we obtain (q) = 0.78. iDe Filippis et all 
(|2005|) also provided the mean axis ratio (q) = 0.81 for 
25 clusters by Chandra and XMM-Newton. These two 
results are consistent with ours ((q) = 0.791 ± 0.004). 

In this paper, we compute the theoretical prediction 



of the PDF based on the phenomcnological model of 
JS02 under the hy drostatic equ i libriu m and isothermal 
assumptions, while. IFlores et aLl(|2007D derived the mean 
cllipticity and dispersion based on 46 haloes in the N- 
body simulation. Their prediction of the mean ellip- 
ticity is (e) = 0.323 ± 0.013 and s £ = 0.138 ± 0.008. 
In order to compare them with our theoretical predic- 
tion, we compute the theoretical value for the same cos- 
mological parameter of their N-body simulation (fi m = 
0.3, A = 0.7, ag = 0.9, h = 0.7). because they used 46 
haloes with M = (1 - 4) x 10 14 h~ 1 M Q , we choose two 
cases of halo mass at z = 0, Mi = 1O 14 /i _1 M0 and 
Mi = 4 x 1O 14 /i _1 M . Performing the Monte Carlo sim- 
ulation (n = 10 5 ) with a = 1 and x'„ = 0.3, we obtain 
(e) = 0.307 and s e = 0.128 for M = M x or (e) = 0.320 
and s e = 0.130 for M — M4. Therefore our model is 
consistent with their model. 

In short summary, our results of the mean axis ratio 
and its dispersion for the observation and the theoretical 
prediction are basically consistent with previous stud- 
ies except for the double peak feature seen in the result 
by Einstein. More detailed comparison needs to reduce 
uncertainty due to the fitting methods or the selection 
effects. 

4. CONCLUSION AND SUMMARY 

In this paper, we have examined the axis ratio distri- 
bution of X-ray halos using the XMM-Newton cluster 
catalogue. By fitting the surface brightness contours by 
ellipses, we confirmed that the radial dependence of the 
axis ratio and the axis direction is relatively small, that 
is, the X-ray halo is well approximated by the triaxial el- 
lipsoid. Constructing the simple model based on the hy- 
drostatic equilibrium of the underlying tr iaxial dark mat- 
ter halo proposed by Uing fc Sutol (|2002| ) , we found that 
the observed PDF of the axis ratio of X-ray halo agrees 
well with the theoretical prediction. We have shown that 
the axis ratio distribution of a sample of X-ray clusters 
predicted by the A CDM model agrees with that mea- 
sured by XMM-Newton, consistent with the findings of 
previous studies using different data sets, which is en- 
couraging given the simplifying assumptions adopted in 
all such model comparisons, including our own. 
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Nishimichi, Thierry Sousbie, and Yasushi Suto for use- 
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(Japan Society for Promotion of Science) Grant-in-Aid 
for science fellows. This work is also supported by Grant- 
in-Aid for Scientific research from JSPS and from the 
Japanese Ministry of Education, Culture, Sports, Sci- 
ence and Technology (Nos. 20-10466 and 22-5467), and 
by the JSPS Core-to-Core Program "International Re- 
search Network for Dark Energy" . 
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Table 1 

Axis ratio q(Ri) and its standard deviation a q (Ri) for the threshold of the contours Ri. 



Z - 1 1 - - ,_.J_ 

Cluster 




Uq(Kl) 


g(R2) 




q(R3) 


o- q (K3) 




<J q (K4 ) 


2A0335+096 


0.84 


0.004 


0.84 


0.01 


0.79 


0.02 






A13 


0.7 


0.06 


0.76 


0.05 


0.75 


0.07 






A68 


0.69 


0.1 


0.67 


0.04 


0.73 


0.05 






A85 


0.94 


0.02 


0.83 


0.02 










A133 


0.86 


0.03 


0.84 


0.02 


0.8 


0.05 






A209 


0.65 


0.06 


0.8 


0.05 


0.79 


0.05 


0.79 


0.09 


A262 


0.8 


0.02 














A383 


0.91 


0.03 


0.93 


0.03 










A399 


0.89 


0.04 


0.91 


0.04 


0.69 


0.03 






A400 


0.67 


0.05 














A478 


0.74 


0.005 


0.77 


0.01 


0.84 


0.01 


0.88 


0.03 


A496 


0.82 


0.01 


0.77 


0.02 










A576 


0.81 


0.04 


0.89 


0.03 










A665 


0.73 


0.04 


0.82 


0.04 


0.94 


0.03 


0.96 


0.04 


A773 


0.72 


0.04 


0.88 


0.05 


0.77 


0.05 






A1068 


0.7 


0.03 


0.72 


0.03 


0.75 


0.06 






A1413 


0.74 


0.03 


0.68 


0.02 


0.74 


0.03 


0.78 


0.04 


A1589 


0.66 


0.11 


0.67 


0.06 


0.62 


0.05 






A1650 


0.73 


0.01 


0.78 


0.02 


0.82 


0.03 


0.73 


0.05 


A1689 


0.9 


0.03 


0.87 


0.02 


0.84 


0.04 


0.82 


0.04 


A1775 


0.81 


0.02 


0.94 


0.02 


0.87 


0.07 






A1795 


0.76 


0.004 


0.79 


0.01 


0.76 


0.01 






A1835 


0.94 


0.02 


0.9 


0.02 


0.91 


0.03 






A1837 


0.7 


0.01 


0.64 


0.04 










A1914 


0.75 


0.02 


0.85 


0.03 


0.88 


0.04 






A 1 001 


U.o 


U.Uo 


U.o 


U.Uo 










A2029 


0.76 


0.01 


0.83 


0.02 


0.84 


0.03 






A2052 


0.79 


0.01 


0.83 


0.01 










A2065 


0.8 


0.02 


0.71 


0.02 


0.64 


0.02 


0.63 


0.04 


A2163 


0.65 


0.06 


0.88 


0.04 


0.74 


0.03 


0.72 


0.03 


A2199 


0.81 


0.01 














A2204 


0.9 


0.01 


0.87 


0.02 


0.83 


0.03 






A2218 


0.8 


0.05 


0.85 


0.03 


0.81 


0.04 






A2319 


0.7 


0.01 


0.77 


0.02 










A2589 


0.73 


0.01 


0.69 


0.03 
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Table 2 

continued. 





g(Ri) 




g(R 2 ) 


<?a(Ra) 


q(R 3 ) 


Va(R'\) 


l(Ri) 




A2597 


0.82 


0.01 


0.79 


0.01 


0.85 


0.02 


0.7 


0.04 


A2626 


0.83 


0.03 


0.88 


0.02 


0.77 


0.04 






A2667 


0.73 


0.03 


0.83 


0.03 


0.81 


0.05 






A2717 


0.98 


0.02 


1.0 


0.02 










A3112 


0.76 


0.01 


0.79 


0.02 


0.76 


0.03 






A3158 


0.72 


0.03 


0.79 


0.02 


0.67 


0.02 






A3558 


0.67 


0.01 


0.77 


0.01 










A3560 


0.91 


0.05 


0.92 


0.05 










A3581 


0.81 


0.01 


0.79 


0.07 










A3827 


0.9 


0.02 


0.93 


0.03 


0.91 


0.03 






A3911 


0.27 


0.12 


0.59 


0.02 


0.64 


0.02 






A3921 


0.73 


0.02 


0.73 


0.02 


0.67 


0.02 


0.59 


0.02 


A4059 


0.76 


0.01 


0.74 


0.02 










AWM7 


0.69 


0.01 














E1455+2232 


0.9 


0.02 


0.9 


0.04 


0.88 


0.05 






EXO0422 


0.9 


0.01 


0.88 


0.03 










Hydra 


0.94 


0.03 


0.89 


0.03 


0.78 


0.03 


0.78 


0.04 


Klcmola44 


0.71 


0.01 


0.66 


0.01 










MKW3S 


0.7 


0.01 


0.74 


0.02 


0.73 


0.02 


- 


- 


MKW4 


0.83 


0.05 














PKS0745-19 


0.76 


0.01 


0.83 


0.02 


0.77 


0.03 






RXCJ0605.8-3518 


0.88 


0.03 


0.72 


0.03 


0.86 


0.05 






RXJ1347-1145 


0.73 


0.02 


0.84 


0.03 


0.84 


0.04 


0.85 


0.05 


Scrsic 159-3 


0.82 


0.01 


0.83 


0.01 


0.83 


0.03 






Triangulum 


0.83 


0.01 


0.7 


0.02 










ZW3146 


0.81 


0.01 


0.85 


0.01 


0.98 


0.03 


0.89 


0.04 




<}(Ri) 




g(R2) 




g(R 3 ) 


<rrr(R 3 ) 


g(R 4 ) 


o-,(i? 4 ) 




0.78 


0.005 


0.81 


0.004 


0.79 


0.006 


0.78 


0.013 



